Global dynamics of cosmological scalar fields - Part I 
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Abstract 



We investigate the Liouvillian integrability of Hamiltonian systems describing a universe 
filled with a scalar field (possibly complex). The tool used is the differential Galois group 
approach, as introduced by Morales-Ruiz and Ramis. The main result is that the generic 
systems are non-integrable, although there still exist some values of parameters for which 
integrability remains undecided. We also draw a connection with chaos present in such 
cosmological models. The first part of the article deals with minimally coupled fields, and 
the second treats the conformal couping. 

1 Introduction 

Homogeneous and isotropic cosmological models, although very simple, explain the recent 
observational data very well [27, 25J. Their foundation is the Friedmann-Robertson- Walker 
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(FRW) universe, described by the metric 



i 2 2 

as = a 



:1) 



where a is the scale factor, d 2 f^2 is the line element on a two-sphere, and we chose to use the 
conformal time rj. As the scale factor depends only on time, so do all the quantities describing 
the matter filling such a universe. These include the density, pressure, a scalar or vector field, 
and of course the cosmological constant (with a rather trivial dependence on rj). 

That last element provides an explanation for the current accelerating expansion of the 
universe (A Cold Dark Matter model [7]), but a better solution still is sought for. A scalar 
field which could model a matter with negative pressure is extensively used for that purpose. A 
Bose-Einstein condensate of a possibly axionic field is one idea [9] , a phantom field violating the 
energy principle is another [1J. Finally, scalar field could also be the mechanism behind inflation 

A universe filled with only one component seems simple enough but it is not the case here. 
Chaos has been studied in conformally coupled fields by means of Lyapunov exponents, pertur- 
bative approach, breaking up of the KAM tori [5J [3], non-integrability with Painleve property 
Chaotic scattering has been found in minimally coupled fields [SJ, as has dynamics [20J. 

The Hamiltonian of these systems has indefinite kinetic energy part, and to cast it into 
a positive-definite form a transition into imaginary variables is used. It has been done for a 
conformally coupled field |12| . but some authors argue that [21] there are physical limitations 
which forbid extending the solutions through singularities such as a = 0, and an imaginary scale 
factor seems even less realistic. 

Despite the fact, that systems of the considered type were called non-integrable, there was 
no rigorous proof of that proposition. However, the Liouvillian integrability can be studied 
successfully, as we try to show in these two articles. 



2 Liouvillian integrability 

Here we only describe the basic notions and facts concerning the Morales-Ramis approach 
following [214 l30l I3"T] and differential Galois theory following [21 [2S]. A more detailed and 
complete presentation of Morales-Ramis theory can be found in [6j I2T1 [23] and of differential 
Galois theory in [T3l [2D [2S] . 

Let us consider a system of differential equations 

x = v(x), teC, x £ M, (2) 

defined on a complex n-dimensional manifold M. If <p(t) is a non-equilibrium solution of ([2j), 
then the maximal analytic continuation of <p(t) defines a Riemann surface T with t as a local 
coordinate. Together with system ([2]) we can also consider linear variational equations (VEs) 
restricted to TrM, i.e. 

i = T(v)£, T(v) = ^p^, £eT r M. (3) 
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We can always reduce the order of this system by one considering induced (also linear) system 
on the normal bundle N := T^M/TT of T [16] 



7? = ^(T(^)(^- 1 0) 



rj £ N. 



(4) 



Here it : TpM — > N is the projection. Obtained in this way, a system of s = n — 1 equations is 
called the normal variational equations (NVEs). If system ([2]) is Hamiltonian then n = 2m and 
the order of VEs can be reduced by two. Namely using a first integral - the Hamiltonian we 
can restrict system ([2]) to the level H(x) = E. Next the obtained system of n — 1 differential 
equations is reduced once more by one using the described above reduction on the normal bundle. 
As result we obtain NVEs of order s = 2(m — 1) which are additionally Hamiltonian. 

An analytic continuation of S(i) - a matrix of fundamental solutions - along a loop 7 gives rise 
a new matrix of fundamental solutions S(t) which not necessarily coincides with S(i). However, 
solutions of a linear system form n dimensional linear space, so we have 3(t) = E(t)M 7 , for 
certain nonsingular matrix M 7 E GL(n, C) which is called the monodromy matrix. Taking into 
account loops in different homotopy classes, this leads to the monodromy group Ai. 

The monodromy group of NVEs was used by Ziglin to formulation necessary conditions for 
the existence of maximal number of first integrals (without involutivity property) for analytic 
Hamiltonian systems [TBI ED., 31J . The Ziglin method was applied to an analysis a number 
various system: e.g the rigid body [31], the Henon-Heiles system [31] etc. However its range of 
possible applications is strongly restricted by difficulties in calculation the monodromy groups. 
These groups are known for a very limited number of equations. 

Recently Morales-Ruiz and Ramis generalised the Ziglin results replacing the monodromy 
group M. by the differential Galois group Q of NVEs. There are two important advantages of 
this replacement 

• always A4 C Q and integrability theorems formulated by means of Q are "stronger", 

• Q has better algebraic properties and there are tools to calculations Q for some classes of 
differential equation. 

The main problem of differential Galois theory is a problem of solvability of linear differential 
equations and systems of linear differential equations in a class of "known" functions. 

Let F be a differential field of characteristics zero, i.e. a field with an additive operation 
' : F — > F called derivation satisfying the Leibnitz rule. In this paper we meet two examples: 
F = C(z) - the set of rational functions with coefficients in C and M(z) - the set of meromorphic 
functions on Riemann surface T. In both examples the derivation is the standard differentiation 
d/dz. The kernel of derivation is a sub-field of F called the sub- field of constants and we denote 
it as F'. In both our examples F' = C. 

Let F and L are differential fields. The field L is a differential field extension of F if the 
derivation on L restricted to F coincides with that of F and we denote F C L. 

What we consider is a system of linear equations 



X' = AX. 



d_ 

dz 



(5) 
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where A is a matrix of dimension s with entries in F and X = [x±, . . . ,x s ] . This system can 
be always transformed to the form 

yM + a s - iy ( s -V + ... + ai y' + a y = Oi, a { G F. (6) 

There is also a direct formulation of differential Galois theory for systems of linear differential 
equations but a majority of literature concerns a linear differential equation of degree s. For 
this reason we transformed system §J)§ into ©• 

Let y^i), ■ ■ ■ ,U(s) are linearly independent over F' solutions of ©. Usually yi G" F and we 
construct the smallest differential extension K containing all solutions of (jBJ) and such that 
K' = F' . It is called a Picard-Vessiot extension of F. We can imagine K as a set F with 
variables z, y{\), • • • , V( s ), y'n\, ■ ■ ■ , y'r a y ■ ■ ■ > Vnl, ■ ■ ■ , y[*j- If F has the characteristic zero and F' 
is algebraically closed, then a Picard-Vessiot extension exist and is unique up to differential 
isomorphism. 

Next we precise what we understand by "known" functions. A solution y of ([6j) is 

1. algebraic over F if y satisfies a polynomial equation with coefficients in F, 

2. primitive over F if y' £ F, i.e., if y = fa, for certain a £ F, 

3. exponential over i* 1 if y'/y G -F, i.e., if y = exp J a, for certain a £ F. 

We say that a differential field L is a Liouvillian extension of F if it can be obtain by successive 
extensions 

F = K C Kx C • • • C K m = L, 

such that Ki = i^_i(yj) with y^ either algebraic, primitive, or exponential over Ki—\. Our 
vague notion "known" functions means Liouvillian function. We say that ([6]) is solvable if its 
Picard-Vessiot extension is a Liouvillian extension. 

Using this definition we cannot check solvability of ([6|) in a class of Liouvillian functions 
because in general we do not know explicit forms of solution. But solvability of © translates 
into properties of the differential Galois group of Q. This group can be defined as follows. For 
Picard-Vessiot extension K D F we consider all automorphisms of K (i.e. invertible transforma- 
tions of K preserving field operations) which commute with differentiation. An automorphism 
g : K — > K commute with differentiation if g(a') = (g(a))' , for all a G K. The set of all such 
automorphisms we denote A. Let us note that automorphisms A form a group. The differential 
Galois group (DGG) Q of extension K D F, is, by definition, a subgroup of A such that it con- 
tains all automorphisms g which do not change elements of F, i.e., for g G Q we have g(a) = a 
for all a G F. 

Equivalent linear systems of differential equations have the same differential Galois groups 
|26j and differential Galois group of system ([5]) is isomorphic to a differential Galois group of 
associated equation © of order s. 

Let us list basic facts about the differential Galois group 

• Q is an algebraic subgroup of GL(s,-F'). Thus, it is a union of disjoint connected compo- 
nents. One of them containing the identity is called the identity component of Q and is 
denoted by Q°. 
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• M C Q and in a case when ([6]) is Fuchsian, then A4 is dense in Zariski topology in Q. 

• Every solution of equation is Liouvillian iff Q° conjugates to a subgroup of triangular 
group. This is the Lie-Kolchin theorem. 

Morales-Ruiz and Ramis showed a connection between integrability of Hamiltonian systems 
and properties of the identity component of the differential Galois group of NVEs |2H 123] . 

Theorem 1 Assume that a Hamiltonian system is meromorphically integrable in the Liouville 
sense in a neighbourhood of a particular solution. Then the identity component of the differential 
Galois group of NVEs is Abelian. 

Application of the Morales-Ramis theory consists of 

• finding of particular solution(s), 

• writing of variational equations and separating normal variational equations, 

• checking of the identity component of differential Galois group of normal variational equa- 
tions. 

In applications, the last step is the most difficult and successful work depends on the form of 
NVEs. 

Sometimes the application of variational equations to integrability studies is insufficient. 
Then, often, variational equations of higher orders are useful. For their definition and applica- 
tions in integrability theory consult for example |18j . 



3 Minimally coupled field's setup 

The standard, general action for a complex field ifi is given as 



-gd x, 



(7) 



where 1Z is the Ricci scalar, A the cosmological constant, and m the so called mass of the field. 
In the FRW model the metric depends only on time, and so must the field. We choose to use 
the conformal time rj, so that the Lagrangian becomes 



C = 6(a"a + Ka 2 ) + ^\iP'\ 2 a 2 



m 



2Aa 4 , 



(8) 



with the prime denoting the derivative with respect to time, and K the index of curvature. We 
also dropped a coefficient which includes some physical constants and the part of the action 
related to the spatial integration. 

Next we subtract a full derivative 6(a'a)', and use the polar parametrisation for the scalar 
field tp = a/12 0exp(z#) to get 



C = 6(Ka 2 - a' 2 ) + 6a 2 (0 /2 + tfO' 2 ) - 6^UV 



2Aa 4 



(9) 
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and obtain the Hamiltonian 



H = k ~ pl + a^ P °) ~ 6Ka2 + 2A ° 4 + V a ^ 2 - ^ 

Since 6 is a cyclical variable, the corresponding momentum is conserved so we substitute p 2 , = 
2u 2 . To make all the quantities dimensionless, we make the following rescalings 



m 2 



m 2 h 2 \K\, A^3L\K\, uo 2 - 72u; 2 |K|, p 2 a ^ 72p 2 a \K\, p 2 ^ 72p 2 \K\. (11) 



There is no need of changing the variables a and 4> along with their momenta, as this is really 
changing the time variable rj, and thus the derivatives to which the momenta are proportional. 
This result also in dividing the whole Hamiltonian by 6y2|if | to yield 

V2H = \ {-p\ + -U) - ^a 2 + La 4 + r»W + ^. (12) 

If the spatial curvature is zero, any of the other dimensional constants can be used for this 
purpose, so without the loss of generality we take the right-hand side to be the new Hamiltonian 

H = \ {-Pi + ^U) - ka 2 + La 4 + m W + (13) 

and in all physical cases, k G {—1, 0, 1}, cu 2 ^ 0, m 2 S 1, L S 1, = 0. We extend the analysis 
somewhat assuming that the Hamiltonian might be equal to some non-zero constant £sR. We 
will later see, that our analysis includes also the possibility of these coefficients being complex. 

Note that for a massless field, the system is already solvable, as shown in appendix A. 

From this point on, we take uj = 0, which means the phase is constant. Since the model has 
U(l) symmetry, we can always make such a field real with a rotation in the complex ip plane. 

Under this assumption the Hamilton's equations of system (fT3l) are 



-Pa, Pa = 2ka - 4a 3 (L + m 2 (p 2 ) + — p\ , 



a 



1 

° 3 (ID 

2 P<Pi P<t> = -^m^a 1 



We note that there is an obvious particular solution, which describes an empty universe: 
4> = P(f> = 0, a = q, p a = —q. Thanks to the energy integral E = \q 2 + kq 2 — Lq 4 , it can be 
identified with an appropriate elliptic function. 

4 Non-integrability in the A = case 

The system now has the following form 



-Pa, Pa = 2ka - Am 2 a 3 (j) 2 + -»pi, 



a3 J 



~2P<t>' P4> = - 2m - 2 a 4< i 



(15) 
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Using the aforementioned particular solution, for which the constant energy condition becomes 
h 



E=\q 2 + kq 2 , we have as the variational equations 



Pa 

V if I 



( 





-1 








\ 




2k 
























i 

? 




V 








— 2m 2 q 4 





/ 



Pa 



(16) 



The normal part of the above system, after eliminating the momentum variation p^P , and writing 
x for 0W, is 



qx + 2qx + 2m 2 q 3 x = 0, 



(17) 



which we further simplify like before by taking z = q as the new independent variable, and using 
the energy condition to get 



z(E - kz 2 )x" + (2E - 3kz z )x' + mVx = 



2\_/ 



,2„3 „ 



(18) 



We check the physical E = hypersurface. This requires k ^ for otherwise the special solution 
would become an equilibrium point. Introducing a new pair of variables 



we finally get 



w(s) 



d 2 w 



w 



2 m 



x(z), 



ds 2 



1 K 4/i 2 - 



1, 



)w, 



(19) 



(20) 



with fj, = ±1, and k = 0. This is the Whittaker equation, and its solutions are Liouvillian if, 
and only if, (k + \i — |, k — fj, — 5) are integers, one of them being positive and the other negative 
[21j . As this is not the case here, this finishes the proof for k ^ 0. 

Non-integrability on one energy hypersurface means no global integrability, for the existence 
of another integral for all values of E would imply its existence on E = 0. However, there might 
exist additional integrals for only some, special values of the energy. It is straightforward to 
check with the use of Kovacic's algorithm [15] . that this is not true here. In cases 1 and 2, there 
is no appropriate integer degree of polynomial needed for the solution, and case 3 cannot hold, 
because of the orders of the singular points of the equation. 

If k = 0, a change of the dependent variable to w(z) = zx(z), reduces equation (|18p to 



Ew" + m 2 z 2 w = 



(21) 



which is known not to posses Liouvillian solutions |15| . 

We notice that when L = E = k = 0, the system can be reduced to a two-dimensional one. 
In fact, the reduction is still possible when L ^ 0, so we choose to present in the next section. 
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5 Non-integrability in the A ^ case 

We use the nonzero constant L to rescale the system as follows 



9i 

92, 



Pa 



Pi 
P2 



(22) 



so that the equations become 

Qi = - 



Pi, 



1 

92 = — oP2, 

91 



p 2 = -2bq 2 q 1 , 



9i 



(23) 



where b = m 2 /L. The energy integral, for the previously defined particular solution, now reads 
£ = EL = \q 2 + kq 2 — q A , where q has been rescaled according to (|22|) . 
As before, we are interested in the variational equations, which read 



9« 

V p 2 1} / 



/ 

2(fc - 6g 2 ) 


V o 



-1 








\ 


/ 


9r 1} 















P? 










i 

¥ 






$ 







-2bq 4 





/ 


V 




) 



(24) 



and writing x for , and y for p$p , the normal part is 

1 

x = —y, 



Or alternatively 



-2bq A x. 



x + 2-x + 26g x 
9 



(25) 



(26) 



5.1 £ = 

We first pick the particular solution of zero energy, as global integrability implies integrability 
for this particular value of the Hamiltonian. It is important to remember, however, that the 
converse is not true. 

The normal variational equation is cast into rational form by changing the independent 
variable to z = q 2 /k (for k ^ which implies k 2 = 1), and using the energy integral 



a 5z- A . 
x + — —x + 



2z(z - 1) Az{z - 1) 



x = 0, 



(27) 
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with the respective characteristic exponents 



z = 0, p = -l,0 

z = ^ p = ^\ (28) 
z = oo, p = -(3 - V9 - 46), -(3 + V9 - 46) 

By Kimura's theorem [J3], the solutions of equation (|27j) are Liouvillian if, and only if 9 — 46 = 
[2p — l) 2 , p £ Z. As before, this means that for global integrability this condition must be 
satisfied. 

For k = the solution of NVE is x\i = q~ 2prxl > 2 , and the reduction to two degrees of freedom 
is possible, as mentioned before. 

5.2 E^O 

The special solution, is now directly connected to the Weierstrass p function, for if we 
introduce a new dependent variable v with 

e = \v + ^ (29) 
the energy integral implies that it satisfies the equation 

v 2 = 4v 3 - g 2 v-g 3 , (30) 

where 

92 = y (k 2 - 3£), g 3 = ^k(2k 2 - 9S), (31) 

and the discriminant A = W24£ 2 (k 2 — 4£), which we take as non-zero to consider the generic 
case, 
reads 



case. Thus, taking w = q^ q, and eliminating p% as before, the normal variational equation 



w = [Ap(t;g 2 ,g 3 ) + B]w, (32) 

with A = 2 — b and B = — |fc(l + 6). This is the Lame differential equation, whose Liouvil- 
lian solutions are known to fall into three mutually exclusive cases, which are exactly those of 
Kovacic's algorithm: 

1. The Lame-Hermite case, with A = n{n + 1) = 2 — 6, n £ Z+. This implies that 9 — 46 = 
(2n + l) 2 . Since 6 does not limit the values of n in any other way, we have to use the 
higher variational equations, to prove that the Galois group is non-abelian. To do that, it 
is convenient to change the variables of equations ([23D in the following way 



qi = wx, pi = -w 2 , 

w 3 (33) 
q 2 = —, p 2 = W!W 4 - w 2 w 3 , 
wl 
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so that the variational equations are all Lame's equations 



;,(!) 



w 



IV. 



w 



;,0) 



V 



;,(!) 



w 



1 \ 

Aip(7/) + Bi 

1 

A 2 p(v) + B 2 J 



where p{rj) is the one given by equation ([50]) . and 



w. 



w 



w 



(1) 

2 

(1) 
3 

(1) 



(34) 



M =6, 



2k 



A 2 =n(n + l), B 2 = -k(n 2 + n-3). 



(35) 



Using the same method as in [18] . we find that the fourth order variational equation's 
solution involves a logarithm, provided that k ^ 0, thus rendering the Lame-Hermite case 
non-integrable. Note, that HVE can only be used in the first case of Lame's equation [21] . 

2. The Brioschi-Halphen-Crawford case, where necessarily n is half an integer, i.e. n + \ = 
I G N, and as before 9 - 46 = (2n + l) 2 = (2/) 2 . 

3. The Baldassarri case, with n+ \ S ^> an d additional algebraic restrictions 
on B,g 2 , and 53. However, inspecting the third case using Kovacic's algorithm, we see 
that 9 - 46 = (2ra + l) 2 must be a square of an integer which is impossible if n belongs to 
the family mentioned. 



5.3 E = k = 

As mentioned in the A = section, we can transform the system to a two-dimensional 
one. In order to do that, time needs to be changed from the conformal to the cosmological one 
dr/ — * dt = adr], in the original equations ([14]) . We then take as the new momenta the Hubble's 
function and the derivative of <p 

h — -— - 

a dt a 2 /o^\ 
_ d(j) _ 

W '~ dt ~ a 3 ' 

(This to is not to be confused with the one introduced in section 3.) Accordingly we have 

da 

— = ah, 



dt 

# 

dt 
dh 

dtv 



OJ, 



(37) 

4L + im 2 (p 2 -uj 2 -2h 2 , 
—2m 2 d> — 3ujh. 
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Thus, we are left with a dynamical system in the (h, (ft, uj) space, as a decouples. Furthermore, 
the energy integral is now 

= ia 4 (2L + w 2 + 2mV -/i 2 ), (38) 

so for a(t) which is not trivially zero, it gives a first integral on the reduced space. Choosing an 
appropriate variable a, suggested by the form of this integral 



Vh 2 - 2L . . 

= sm a , 

V2m V ; ' (39) 



uj = \/h 2 — 2L cos(a), 

we finally obtain 

— = v2m + 3/isin(a) cos (a), 

t (40) 

H£ = -3(/l 2 -2L)cOs(a) 2 . 

The problem of such reduction was also discussed in [10]. It is argued that there can be no 
chaos in this system, but its integrability - which would be one more first integral - remains 
unresolved. 



6 Conclusions 

The main result of our paper can be summarised as follows. 

When A = 0, and if the system is meromorphically integrable, then necessarily E = k = 0, 
in which case it is possible to reduce it to two-dimensions, and there is no chaos. 

When A 0, we consider the physical E = hypersurface first. If the system is integrable, 
then either k = (reducibility as above), or 9 — 4m 2 /L = (2n + l) 2 , n£Z, 

When A ^ 0, and we consider a generic energy hypersurface, integrability implies that either 
9 — Am? /L = (2n + l) 2 and k = 0, or 9 — 4m 2 /L = (2n) 2 (regardless of the value of k), for some 
n G Z. 

These are, however, only necessary and not sufficient conditions, so that the system might 
still prove not to be integrable at all. In particular, the numerical search for chaos suggests both 
the lack of global first integrals, and crucial differences in the behaviour of the system for real 
and imaginary values of the variables. This might be a clue, that the system might have first 
integrals which are not analytic, and thus not prolongable to the complex domain. A system 
with similar property was studies by the authors in |19| . 

It might also prove useful, to investigate the existence of asymptotic first integrals, but this 
is more in connection with scattering dynamics of conformally coupled fields, which we study in 
the second part of the paper. 
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Appendix A. A solvable case 

When m = 0, the Hamilton-Jacobi equation for the main Hamiltonian fjl3[) will become 
separable, because it can be written as 

„ 2 1 fOW\ 2 1 2 (dW\ 2 , 4 T 6 u 2 

with the full generating function S 1 = W — -Ery. Assuming W = A{a) + F(<p), equation (|4ip can 
be solved with 



A{a) = J J 2 (^La 4 - ka 2 - E + ^ J da. 



The first equation of motion can then be deduced from 



<9£ 

which can be rewritten as 



Or, introducing a new variable v = a 2 , as 

2 

v dr/ 

so that the general solution is 

1 . A; 



where 

52 = y k 2 + 16L£ 

«73 = f^ + gfe 3 -32L 2 J, 



(42) 



dW f da 

rj = I — ; = const, (43) 

2(La 4 -ka 2 - E + ■£) 



f^Y = 2 (La 4 -ka 2 -E + ^). (44) 



^) = 8(Lv 3 — kv 2 — Ev + J), (45) 



a 2 = v = 7^-p(??-%;32,33) + ^-, (46) 



(47) 



and r/o is the constant of integration. Of course, for L = 0, equation (|45j) admits solutions in 
terms of circular functions. 
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The equation for (p(rj) is the following 



dW _ f dcj) f da 

dJ 



/da , , . 
= = const, (48) 
nl. Of T.ni — hr,2 _ TT. _l_ 2. \ 



2(J~f) J a^2(La*-ka2-E + ^) 
which we simplify by using the just obtained solution for v {rj) to get 



const = VJt 2 -"\ f*n (49) 

As v is an elliptic function of order two, the second integral can be evaluated by means of the 
Weierstrass zeta function to yield 



const = V2J + WTj [Cim) " CM] v + 111 



(50) 



where 771^, are the zeroes of v(rj), given by 

3p(77i, 2 ;S2,53) = -2k, (51) 

and the constant of integration can be determined from the boundary conditions on the field (j). 
Again, for J = 0, the integrals in (f4*8l) reduce to simpler functions. 
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